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Abstract 

We derive a closed-form, analytical expression for the spectrum of long- 
wavelength density perturbations in inflationary models with two (or more) 
inflaton degrees of freedom that is valid in the slow-roll approximation. We 
illustrate several classes of potentials for which this expression reduces to a 
simple, algebraic expression. 
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I. INTRODUCTION 



Inflationary cosmology, in addition to resolving the cosmological flatness and hori- 
zon problems of the standard big bang model,0 predicts a remnant spectrum of density 
perturbations.ihi The perturbations may have seeded large-scale structure formation and 
may have left an imprint on the cosmic microwave background (CMB) anisotropy. Hereto- 
fore, considerable attention has been given to precise derivations of the perturbation spec- 
trum when a single scalar (inflaton) field drives inflation and simultaneously this field is 
responsible for origin of initial inhomogeneities.i^'ifl Recently, we have shown how differ- 
ent methods of computing the perturbation are related and compared their accuracy.! These 
results have been applied to obtain predictions of CMB anisotropy and large-scale structure. 

In this paper, we derive the explicit, analytical solution for the perturbations in the 
models with two (or more) inflaton degrees of freedom in a case of general potential for the 
scalar fields. To obtain this expression, we must assume that the slow-roll approximation 
for the evolution of the scalar fields is valid. With two or more inflaton fields, the pertur- 
bation spectrum has additional contributions which do not occur in the single-field case; in 
particular, the single field spectrum consists purely of adiabatic fluctuations whereas the 
multi-field spectrum generically has an entropic contribution as well, as discussed for some 
specific cases previously.iHniilii§@ Previously, a formalism has been developed for cases 
in which the two fields are not coupled or have very simple couplings.0'lii'00 However, in 
more complicated cases, only heuristic arguments based on single-field inflation have been 
applied. Yet, the single-field case is unusual because the evolving expectation value of the 
one field serves as a clock that determines when the universe exits inflation and returns to 
Friedmann-Robertson- Walker expansion. This was the basis, for example, of the time-delay 
formalism introduced by Guth and Pii for computing the perturbation spectrum for single- 
inflaton models. With two or more fields, fluctuations in one field can affect the evolution of 
the other field, and the complex conditions under which inflation ends cannot be expressed 
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in terms of one degree of freedom {e.g., some linear combination of fields). Consequently, 
the heuristic arguments are suspect. The issue has become more important in recent years 
because intriguing new models of inflation have been proposed which entail two or more in- 
flaton degrees of freedom: e.g., extended and hyperextended inflation,0'Ei hybrid inflation,lii 
and supernatural inflation.0'0 

For these reasons, it has become essential to have a formalism that applies to the multi- 
field case. Below we develop a procedure based on the natural generalization for the single- 
field case. We assume the slow-roll approximation in which the inflation kinetic energy is 
negligible compared to its potential energy during inflation. We first review the single-field 
case, then a simple case with two decoupled fields, and finally we solve the equations for per- 
turbations in the general case. The central result is in Section IVb, Eqs. (^6[) through (^9D, 
the closed-form expressions for the perturbation spectrum. We present a series of potential 
forms for which the closed-form expressions reduce to simple algebraic expressions. 

II. PERTURBATIONS IN SINGLE-FIELD INFLATION 

We consider first the case of a single scalar field with potential V{(j)). Assuming the 
spatial part of the energy-momentum tensor is diagonal, the metric in longitudinal gauge is: 



where a is the Robertson- Walker scale factor and $ is the gravitational potential. The 
perturbed Klein-Gordon equation which describes the evolution of perturbations in is (in 
units with h = 4itG = c = 1): 



ds 



,2 




(1) 



S(f) + 3H6(j) - \v^6(f) + V"6(j) - 40$ + 2V'$ = 0, 



(2) 



a 



and the perturbed — i Einstein equation is 



$ + i7$ = 00^0, 



(3) 



3 



where (f) — (po + (^0, dot means d/dt, and prime is used for 9/90. Throughout we use 
dimensionless units where AttG = 1, where G is Newton's constant. 

To find the nondecaying solution for the long- wavelength inhomogeneities (for which the 
spatial derivatives in the equations of motion can be neglected) in slowroU approximation, 
the terms proportional to $ or depending on second derivatives in 6(j) can be dropped, 
resulting in the simplified equations 

3H5(I) + V"5(t) + 2V'^ = (4) 

= ^o6(l). (5) 

The (unperturbed) background equations are: 

3H<j)o = -V 

= IV (6) 

H = -4 

where the slow-roll approximation has been assumed in dropping 0o in the first equation 
and 0^ in the second. 

If we introduce a new variable, x = S(f)/V', the perturbative equations reduce to: 

3Hx = -2$ 

(7) 

Substituting the second expression into the first and using the background equation to 
express H in terms of V^, we obtain 

V 

X = --X, (8) 

whose solution is 

X = ^, (9) 
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where C is an integration constant. Therefore, using the definition of x and Eq. (j^), we 
obtain 

50 = V'x = C^ = -2C^. (10) 

To fix the integration constant, we use the quantum de Sitter fiuctuation result that 6(f)k ~ H 
evaluated at horizon-crossing, k = aH. Solving for C above, we obtain 



c^i-H:^ . (11) 



k=aH 



According to Eq. (^, $ = Vx/H. Since our solution is x = C /V , we obtain 

The last expression is obtained by noting = and H/H = 2V/V where, applying the 
slow-roll approximation, we have ignored the infiaton kinetic energy contribution to total 
energy density during inflation. At the end of inflation, H = 0{1)H'^, which means that 

^f^O{l)C^O{l)(H^-i] , (13) 



\ / k=aH 

the standard, lowest-order result.! The fluctuations described in this relation are adiabatic. 
For the cosmic microwave background (CMB) anisotropy in large angular scales, the gravi- 
tational potential sets the CMB temperature fluctuations, 111 6T/T ^ $/3. 

III. TWO DECOUPLED FIELDS 

Next, consider the case of two decoupled flelds 0i and 02 with potential V^(0i,02) = 
^i(0i) + V2(02)- This case has been considered previously.^ In this case, the Klein-Gordon 
equations for the flelds in the slow-roll, long-wavelength approximation takes the form: 

3H5(j)^ + Vl'd(j)i + 2VI<^ = 

(14) 

3HS(j)2 + V^' 6(1)2 + 2V^^ = 



5 



and 



where V{ = dVi/d(j)i and = dV2/d(f)2. 

The (unperturbed) background equations are: 

SH^i = -VI 

3H^2 = -V^ 

= m+V2) 



a2 I i2\ 



Similar to the single-field case, we introduce two new variables Xi and X2 via 

(501 = V^xi and 502 = V"2'x2. 
Then, the system of perturbed equations, Eq. ([T^|T5|), can be rewritten as 

3Hxi = -2$ 

3Hx2 = -2$ 

^{VlXi + V2X2). 

rl 

Subtracting Eq. (|1^) from Eq. ([18|), we obtain xi — X2 = 0, or 

X2= Xi+ D, 

where D is an integration constant. Substituting this into Eq (120|) , one gets that 

^ = ^[{V, + V2)X, + DV2]. 

Using this expression to replace the right-hand-side of Eq. we then obtain the 

form equation for Xi. 
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Vi + V2 DV2 



(23) 



V1 + V2 V1 + V2 

(where the background equation expressing in terms of Vi and V2 has been used). This 
equation may be rewritten as 



^(Xi(Vi + V2)+ W2) =0, 



(24) 



and then can be integrated to obtain 



Xi 



C-DV2 



(25) 



V1 + V2 ' 

where C is the integration constant. The corresponding resuh for X2 is 

= (26) 

The integration constants C and D should be fixed by assuming that 501^2 ^ H a.t horizon 
crossing during inflation, k = aH. Using the relation between xi^2 and 501^2, we obtain 



C = — 
2 



D 



H 



1 / 501 502 



k=aH 



V2 5 (pi 

V1 + V2 01 ' Vi + V2^ 



+ 



(27) 



k=aH 



Also, substituting Eq. (^) into the expression for $, Eq. ([2^), we derive 

H V1 + V2 



(28) 



The first term above is the same as the total contribution in the single field case (see 
Eq. ([T^)); as in the single-field case, it can be interpreted as the adiabatic contribution to 
the fluctuation spectrum. The second term is a new contribution due to entropic fluctuations 
which arises whenever there are two or more inflaton fields. The entropic contribution arises 
because two or more components of the cosmic fluid are undergoing different fluctuations 
and evolution. Hence, the entropic contribution is a characteristic feature of multicomponent 
inflation models.0 



IV. GENERAL CASE 

We now generalize our method to models with two scalar fields in which the kinetic 

energy terms have non-linear sigma-model form and the potential V (0, ijj) is an arbitrary 
function of (f),ip : 



-I 



^/l (0, V) 0;a0=" + \f2 (0, V') V';aV''" " V (<f>, i^) 



y/^d'^x (29) 



Kinetic terms of this type occur, for example, in supergravity models with non-trivial Kahler 
potential. 

A. Derivation 

Variation of this action with respect to 0, leads to the equations 

(Mn;a - lh<t>M'''' - lf2,^M'''' + v<P-Q (30) 

(Mn;a - ih^M''"" - lf2,^M-''' + V^^0 (31) 

We consider a Priedmann-Robertson- Walker Universe with small perturbations 

ds^ = (1 + 2$) dt^ - (1 - 2$) 5ikdx'dx^ (32) 
and decompose fields into homogeneous and inhomogeneous components 

= (/.o(t) + (5(/.(x,i),V' = V'o + (^V' (33) 
The equation for the background then take the form 

(00 + 3//0o) /i + /i0o - \h,^^l - \h,4>i^l + = (34) 

Here all functions depend only on the background variables. The equation for %l) can be 
obtained from the above one if we make the following substitutions: /i f2', 4> ip. 
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The first linearized equation for perturbations is 



/i,^ (00 + 3i/0o) + \h. 



(35) 



The second equation is obtained by the above mentioned substitution. Because we have 
three unknown functions, we need a third relation. It is convenient to choose the 0-i Einstein 
equation: 



(36) 



The above equations can be simplified and integrated explicitly for the longwave pertur- 
bations in slowroU approximation. In this approximation the equations for the background 
simplifies to 

3///10O = -V^ 
3HMo = -V^ 



H^ = -V 



and the equations for perturbations take the form 



3H5<t> + ( ^ 



3//5^ + ( ^ 
h 



h J ,v> h 



(37) 



(38) 
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Let us define the variables x and y as 



and 

Jl 



A/ 

h 



(39) 



In terms of tfiese variables, the above equations-of-motion for the perturbations become 



3Hx + 



- x) + 2$ = 



3i?y H nr \ (x - y) + 2$ = 



(40) 



These equations are very similar to Eqs. (|T8|) and (p!9| ) for the decoupled case, except that 
there is here the middle term which vanishes if = Vi{(f)) + V2{ip). Nevertheless, we can 
obtain a closed-form solution. Subtracting the two equations, we obtain 



3H (y-x)- i^X^tllM^ + {y,/h)Ay,/h) \ _ 



(41) 



which can be integrated 



y — X = 'J exp 



3^ 



(42) 



where 7 is a constant of integration. Using the — ? Einstein equation, $ can be expressed 
in terms of x 



(43) 



Eqs. (|40| ) and (^21) can then be used to find an integral expression for x: 



X 



1_ 
V 



' jv<i>lh),^{vjh) , / 

2(K.//.) 



:dt 



(44) 



where 



F = V^exp 



{v^/h)Av,^/f2) , (v:^//2).(K0//i) 



+ 



dt 
3H 



(45) 



(K^//i) (v:^//2) 

Similar expressions can be obtained for y. Using the definitions of x and y (Eq. ( PU[ ) and the 
equations for the background (P^, the expressions can be simplified to the final, compact, 
closed forms given in the next section. 
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B. A General, Closed- form Expression 



The following closed-form expressions are a compact, general representation of the per- 
turbations with wavenumber k in two-field models in the slow-roll approximation, the central 
result of this paper. The perturbed variables 50, Sip, and $ are fc-dependent, but the sub- 
script k has been suppressed for simplicity. 

S^^t)^,tX}Ar'L((^)] Fii.^atX}± (46) 



A Jto \ \ /i / / ,^ /i 



where 



(48) 



(In re-expressing F, an overall constant has been removed and absorbed into the definition 
of 7.) Here to is an arbitrary moment of time; we take to to be the moment of horizon 
crossing tk when k = aH for the given mode k. The limits of integration indicate that the 
integration variable is to be evaluated at the time to and t. 

There are three integrations constants, a, (3 and 7. Evaluating these expressions at 
horizon-crossing and using the fact that 6(j) ^ Sip ^ H at horizon-crossing, the integration 
constants a and (3 can be determined in terms of the inflaton potential and derivatives 
during inflation. The third integration constant, 7, must be taken in such a way as to satisfy 
Eq. (^2|) , from where it follows that (3 = a + '-y. (One must make use of the definitions of x 
and y in Eq. ( p9|) and the fact that, at horizon-crossing, the integral on the right-hand-side of 
Eq. (|4^ ) is l/Vk=aH-) The value of 7 corresponds physically to the amplitude of the entropy 
perturbations. The gravitational potential is expressed in terms of 6(j) and 6tp as 



$ = — 
2 



{\nV).S^+(\nV)^Jij (49) 
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The above formulae can also be rewritten in other forms which can be useful for approxi- 
mate evaluation of the integrals. For instance, in the standard case where the kinetic energy 
terms are canonical (/i = /2 = 1), the expressions above reduce to: 

J to + 



5^ = 7(lnV),^ f -^^dF + /3i\nV)^^ (50) 

J to ^,ip ^.6 



V. APPLICATIONS 

In a surprisingly wide range of cases, the closed-form integral expressions above can be 
simplified. Our first two examples are simple cases which have been studied previously in 
the literature. We use these to show how our general expressions are to be used and to test 
that they reproduce known results. We then apply the method to more general and more 
realistic models. 

1. Example 1: /i = /a = 1, V = Vi (<^) + V2 (-0) 

We have already presented a derivation for the most trivial case where the two inflaton 
fields are decoupled, V = Vi (0) + V2{ip) (see Section III). Here we show that the answer 
can be reproduced by our more general formulae. We detail a few steps to aid the reader in 
becoming familiar with applying our general formulae. 

If = Vi{(f)) + ¥2(1^), then the expression for F{t) in Eq. (^) can be simplified. The 
first integrand in the exponent can be reduced to — (InV) and the second integrand 
to — {In V)^^dip. The two can be combined into a total differential, —d{\nV); as a result, 
F reduces to F(t) = V/Vk=Ha, where Vk=Ha is the value of potential taken at the mo- 
ment of horizon crossing to = tfc. In the first term of Eq. (^6|), the integral reduces to 

- jl{V;/Vk=Ha)dtP = {V2/V),=Ha ' (^2/^fc=//a) SOthat 
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7 



h a + 7 — 

Vk=Ha V V \y 



k=Hc 



and 



(51) 



k=Ha 



V 



which agrees with Eqs. (p!7[) and (^) in our earher derivation. Comparing these formulae 
with Eqs. (|l^ and ( ]25|) we see that the integration constants C and D in Eq. (pSD are hnear 
combinations of 7, a and [3 : 



C = a + 7(y 



A;=_ffa \ / k=Ha 



D 



7 



k=Ha 



(52) 



This identification between the three integration constants and the two coefficients, C 
and requires use of the constraint that a and (3 must satisfy, /3 — a = 7 ; see remarks 
above Eq. ( ^9]) in Section IVB. Also note that the integration constants are functions of 
the wavenumber k which can be related to dip and 5ip at the moment of horizon crossing 
k = Ha using the relations (^61), ( ^Tf ) and taking in these formulae t = tk : 



a 



k=Ha 



2 V '0 



k=Ha 



^=2 



(53) 



^ / J k=Ha 

Note that the above expressions for coefficients are valid in general case. 

Taking into account the above relations (E2|) we see that these expressions for the con- 
stants of integration are in agreement with formulae (27). We refer to this as Example 1 
Other examples where the expressions for the solutions can be significantly simplified are: 
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2. Example 2: /i = = 1, V = Vi (0) V2 (V') 



The integrands in the expression for F{t) in Eq. (^Hj) and in the expressions for 5(j) and 
if) are precisely zero. Hence, F{t) = 1 and 

S<P = a^, = (3^. (54) 

Eq. ( ^91 ) then reduces to a closed-form expression for 



'^ = --J-l^ + P^], (55) 



m ai 



greement with results obtained previouslyEl for this kind of potential. 



3. Example 3: h = h (V') , h = f2{4>), V = Vi {4>)V2 W 
The solution is 

J2dfi + a- 



r / 7 ^2,1/; t f Af 1/3 ^2,V 



{hh)k=Ha^'i h Jtu ^2 /2 

where a, (3 and 7 are given by the formulae (^31). In particular, when fi 
exp {-iIj/iIji) , /2 = 1, V2 (V') = exp {-ilj/ip2) we obtain 



{(3 exp (ip/ipi) -'yexp{ilJk=Ha/'ipi)) 



m a; 



5^ = -f3/ij2 (56) 
greement with the recent result obtained by Starobinskii & Yokoyama.0 
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4. Example 4: f^ = f^ = l, V = K+ Vi (0) V2 W 



In this case we have 

where as before the coefficients a, (3 and 7 can be expressed in terms of 5(f) and 5ip at the 
momemt of horizon crossing via ( |5^ ) . 

5. Example 5: fi = f2 = l, V = A + aVi (</>) + PV2 (ip) + Vi (0) ^2 (V-) 

This case can be reduced to the previous one by the redefinition of the potentials Vi and 
V2 (shifting them by constant terms). 

In summary, the central results of this paper are in Section IVb, explicit analytical 
solutions of equations for the energy density perturbations in the general case of two scalar 
fields with an arbitrary potential. To apply these solutions to a particular inflationary model, 
one should find ffist the analytical solutions which describe the behavior of the background 
and then use the obtained solutions in the formulae we have provided. A significant feature 
compared to the single inflaton field case is the existence of two constants of integration, 
each set by conditions at horizon-crossing during infiation. The two coefficients can be 
chosen so as to correspond to adiabatic and entropic perturbations. The latter is absent in 
the single-field case, but here is shown to be a general feature. In a future publication, we 
shall discuss the observational consequences of these results; namely, how observations can 
be used to distinguish single- from multi-field infiationary models. 

We thank G. Huey for suggestions for improving the manuscript. This research was 
supported by the Department of Energy at Penn, DE-FG02-95ER40893 (PJS), and by the 
Tomalla Foundation (VM). 
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